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I. INTRODUCTION 



For the past decade or so, the study of first-order phase transitions in cosmology has been 
the focus of much interest due to their possible relevance to the physics of the early Universe. 
Some well-known examples are inflationary models Jl| , the quark-hadron transition || and, 
more recently, the generation of the cosmological baryon asymmetry in the electroweak phase 
transition 

In a first-order phase transition, the initial metastable phase decays to the stable phase 
by the nucleation of bubbles larger than a critical size. This decay may be triggered by either 
quantum or thermal fluctuations, depending on how the ambient temperature compares to 
the nucleation barrier ||]. Within a cosmological context, the cooling is provided by the 
expansion of the Universe; the long-wavelength modes of the order parameter responsible 
for the symmetry breaking transition are coupled to the "environment", which is assumed 
to be in local thermal equilibrium at some temperature T. (Here, we are mainly concerned 
with "late" transitions, for which the typical fluctuation time-scales are much shorter than 
the expansion rate.) 

Of great relevance to the understanding of the evolution of the phase transition is the 
determination of the bubble nucleation rate per unit volume. This is a well-known problem 
in classical statistical mechanics, with an long history ||. Phenomenological field-theoretic 
treatments were developed by Cahn and Hilliard ||, and by Langer || in the context of a 
time-dependent coarse-grained Ginzburg-Landau model. Classical homogeneous nucleation 
theory within a field theoretic context has been recently shown by numerical experiments 
to successfully predict the nucleation barrier |7J . In the case of zero-temperature quantum 
field theory, the study of metastable vacuum decay was initiated with the work of Voloshin, 
Kobzarev, and Okun M, and was put onto firm theoretical ground by Coleman and Callan 
in the late seventies 0. Finite temperature corrections to the vacuum decay rate were 
first considered by Linde [HJ, who argued that temperature corrections to the nucleation 
rate are obtained recalling that finite temperature field theory (at sufficiently high temper- 
atures) in d + 1-dimensions is equivalent to <i-dimensional euclidean quantum field theory 
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with h substituted by T. Thus, in d + 1 dimensions, the nucleation rate is proportional 
to exp [—Sg (ipb) /T], where S^((pb) is the <i-dimensional euclidean action evaluated at its 
extremum (specifically, a saddle point), the critical bubble, or bounce, ip b {f)- The usual 
expression for the nucleation rate per unit volume used in the literature is |H| 

S|Mr,T),T)" 



v _ T ( SUMr,T),T) \^ j det[-V 2 + V^ f ,T)} 1' 

I 2ttT ) 1 det'[-V 2 + Mr, T),T)]J P 



where ipf is the value of the field <p at the metastable minimum, and the prime in the 
determinant in the denominator is a reminder that one should omit the zero and negative 
eigenvalues, associated with the translation symmetry of the bubble and with its instability 
(being a saddle point configuration), respectively V e g(<p, T) is the one-loop approximation 
to the finite temperature effective potential, and V^((p,T) = 9 ^jf'^ \<p=y 

There are three important points here. The first is that in order to estimate the de- 
terminantal prefactor (the "equilibrium" part of the prefactor; there is a dynamic factor 
which can not be obtained by using equilibrium arguments) one usually proceeds by in- 
voking dimensional arguments to approximate it by a term of order Tq (Tc is the critical 
temperature). How good is this approximation? Clearly, in most cases it is impossible to 
evaluate the determinants exactly. But can one obtain a better approximation than the 
simple use of dimensional arguments? The second, and most important, point is that the 
critical bubble configuration ipi,{r,T) used to evaluate the nucleation barrier, denoted above 
by Sfi(tpb(r,T)), was obtained from an effective potential which includes corrections coming 
from scalar loops. Hence the temperature dependence in ip(r,T). We will show here that 
this procedure is not in general justified and is only a good approximation if the corrections 
from scalar loops are negligible. Finally, in the expression for the temperature corrected 
barrier in Eq. ( |1.1| ), S^(ifb(r,T),T), one uses the temperature corrected effective potential, 
V e s(<p,T) as opposed to the tree level potential. Thus, it is claimed that S^(ipb(r,T),T) is 
equivalent to the free energy of the temperature dependent bounce, given by 



S%(<p b (r,T),T) = d 3 x 



i(V^(r,T)) 2 + Kfr(^(r,T),T) 



(1.2) 



As far as we know, apart from the work of Affleck in the context of finite temperature quan- 



tum mechanics [fLl]j , this point has never been properly addressed in the literature. How 
did the temperature corrected potential appear in the exponent? Is the exponentiation of 
the massless modes sufficient? In fact, most of the work done on cosmological phase tran- 
sitions in which temperature effects are important (including some by the present authors) 
simply invokes Linde's results. Given the many applications of finite temperature vacuum 
decay in cosmology, we feel that this important question should not be left unscrutinized. 



This concern has also been expressed in recent works by Csernai and Kapusta and by 
Buchmiiller, Helbig, and Walliser ||13fl . Both works attempted to improve on Linde's results, 
by generalizations of Langer's work. Csernai and Kapusta obtained an expression for the 
dynamical prefactor by using a relativistic hydrodynamic approach, while Buchmiiller et 
al. obtained an approximate expression for the decay rate in scalar electrodynamics (and 
more recently, with Z. Fodor, in the standard electroweak model) by integrating out the 
electromagnetic degrees of freedom from the partition function. However, a more detailed 
analysis of the nucleation barrier and how it compares to the usual result is still lacking. 

In this paper we address the three points raised above. We will be mostly interested in 
the regime in which thermal fluctuations are much larger than quantum fluctuations. This 
way we avoid the question of how to match continuously the two regimes, although we believe 



this to be a very important question [14|. (See also Refs. [ I l.| and [15].) The hope is that 
in most situations of interest the transition will be dominated by one or the other regime. 
By a saddle-point evaluation of the partition function in the case of a self-interacting scalar 
field, it is possible to show that the temperature corrections to the nucleation barrier can be 
interpreted as entropic contributions due to stable vibrational modes on the tree-level bounce 
configuration, (fb{r). In other words, the first corrections to the energy of the critical bubble 
configuration come from temperature induced stable fluctuations on the bubble, which will 
modify its volume and surface energies. We will show that these corrections are given by the 
temperature corrected effective action evaluated at the tree-level bounce <Pb{ r )- I n expression 
( |1 . 1|) , the bounce is obtained from the effective potential which includes scalar loops. The 
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difference between the two nucleation barriers will be important whenever scalar loops are 
not negligible. We will show that they become particularly important within the so-called 
thin- wall limit, that is, in the vicinity of the critical temperature for the transition. This is 
perfectly consistent with the fact that large entropic corrections are expected near the critical 
temperature. We will obtain this result by a perturbative evaluation of the determinantal 
prefactor. In principle, the determinantal prefactor can be evaluated in two ways. Clearly, 
the computation can be done directly if we know the eigenvalues related to a given bubble 
configuration. This method is not very useful in practice, since we in general do not know the 
eigenvalues. (Unless, of course, we obtain them numerically.) Writing down explicitly the 
eingenvalue equations, and using a thin-wall approximation to the bubble configuration, we 
show how the temperature corrections to the nucleation barrier originate from fluctuations 
about the critical bubble configuration. Even though the thin-wall approximation is not very 
useful in realistic situations, there is no reason to believe that thicker wall bubbles will behave 
any differently. (Unless the transition becomes too weak, in which case nucleation of critical 
bubbles may not be the relevant mechanism for the transition ||P7|] .) The second approach we 
use to evaluate the prefactor relies on a perturbative expansion of the determinants. Within 
first-order, it is again possible to show how the prefactor accounts for the temperature 
corrections to the nucleation barrier. 

The paper is organized as follows. In Section 2 we briefly review Langer's formalism for 
obtaining nucleation rates, adapted to field theory at finite temperatures. That is, we obtain 
the partition function for the metastable phase plus a nucleating fluctuation by a saddle-point 
evaluation of the functional integral. In Section 3 we show how the determinantal prefactor 
can account for the finite temperature corrections to the nucleation barrier. For simplicity, 
the calculation is performed in the context of the thin-wall approximation for the bubble 
profile, although in principle one could obtain results for any configuration. In Sections 2 
and 3, for the sake of clarity, the discussion is somewhat oversimplified. We assume that the 
system is initially in a metastable state and study only the scalar degrees of freedom in the 
problem. This situation is not unrealistic, as it can be reproduced in numerical simulations 
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of vacuum decay |7j. In Section 4 we study a model of a scalar field coupled to fermions. 
This example is particularly interesting as it illustrates how a thermal state evolves into a 
metastable state due to radiative corrections, very much like in the standard electroweak 
model. (Recall that in a cosmological context the cooling is provided by the expansion of the 
Universe.) A related problem has been recently studied by E. Weinberg, in the context of 
massless scalar models for which symmetry breaking occurs due to radiative corrections. In 
the regime dominated by quantum fluctuations, Weinberg showed how radiative corrections 



induce a metastable state and how it is possible to evaluate its decay rate ]I8| . We show how 
to obtain an effective partition function for the scalar field by integrating out the fermionic 
modes, and proceed to obtain the decay rate. We then compare our results for the nucleation 
barrier to the results obtained using Eq. (\1 . 1| ) . Conclusions are presented in Section 5 and 
two Appendices are included to clarify a few technical points. 

II. FINITE TEMPERATURE DECAY RATE: GENERAL FORMALISM 

Consider a scalar field model with four-dimensional euclidean action 

S E = J d A x E C E , (2.1) 

where Ce is the euclidean lagrangian density given by 

£ s = ^(<V>) 2 + U(0) , (2.2) 

and the potential in ( |2.2| ) has a metastable minimum at = tpf and a stable minimum at 
= !f t , as shown in Fig. 1. [We will only consider potentials with two minima here.] 

Let us assume that the system is prepared initially in the metastable phase, without 
worrying for the moment about how this is done. (See Section 4.) The metastable phase 
will decay into the stable phase by the nucleation of bubbles larger than a critical size. (For 
a review see, e.g., As is well-known, in order to study the decay of the false vacuum at 
finite temperature we impose the periodic boundary condition (anti-periodic for fermions) 



0(0, x) = (f)((3fi,x), so that the euclidean action becomes [|T| 
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S E [(f>] = J dr I (fx 



(2.3) 



The partition function of the system is given by a functional integral over all possible field 
configurations weighted by their euclidean action, 



D(j)e- bEW . (2.4) 

Following Langer |J, we describe the nucleation of bubbles of the stable phase inside the 
metastable phase under the assumption that a dilute gas approximation for these droplets is 
valid. Unless the transition is weakly first-order, this should be a very good approximation 
to describe the early stages of the transition, when bubble collisions and other complicated 
kinetic effects can be neglected. The critical configuration is an extremum of the euclidean 
action, 

^fl^. = 0, (2.5) 
being thus a solution of the equation of motion, 

(J^ + V 2 )0 = n0), (2-6) 

with boundary conditions, lim T ^± oo 0(r, x) = iff, and lim[j|_> C)0 = iff. Coleman, Glaser, and 
Martin [fL9[| , have shown that the configuration with minimum energy, i.e., the minimum 



of all the maxima, will have 0(4)-symmetry. As argued by Linde [fL0| , for sufficiently high 
temperatures the problem becomes effectively three-dimensional, and the saddle point will 
be given by the 0(3)-symmetric, or static, solution of 

d 2( P + 2 <ty = y,,^ , 2 ? , 

dr 2 r dr 

with boundary conditions, lim r _» +00 (f> = ipj and ^| r =o = 0. Note that the potential in ( |2.7| ) 
is the tree-level potential. For future reference we note that when the false vacuum energy 
density [Ay = V((pf) — V(ipt)] is much smaller than the barrier height [Vh, see Fig. 1], the 
bubble radius R is much larger than the wall thickness AR ~ m _1 , where m is a typical 
mass scale in the problem. In this case, the solution to Eq. ( |2.7|) can be estimated by the 
so-called thin-wall approximation and is given by 



(2.8) 



(p t , < r < R - AR 

fwaii(r), R- AR<r < R + AR 
iff, r>R + AR 

which describes a bubble of radius R of the stable phase <p t embedded in the metastable 
phase ipf. (p wa ii(r) describes the bubble wall separating the two phases. 

The main advantage of Langer's approach is that in a dilute gas of bubbles one can infer 
the thermodynamics of the system from the knowledge of the partition function of a single 
bubble. We can write the partition function for the system with a bubble of the stable 
vacuum inside the metastable vacuum as 



Z = Z(<f f ) + Z(ip b ) 



(2.9) 



where Z(tpf) and Z{ip^) are the partition functions of the system for the vacuum field con- 
figuration iff and for the bubble field configuration respectively |20|. The generalization 
of ( |2.9|) for several bubbles is given by 



Z ~ Z{tp t ) + Z{<p f ) 
~ Z((ff) exp 



Z(<fb 



n 2 



2! [Z(<p f )\ 



+ ... 



Z(<Pf) 



(2.10) 



The proof of ( p.lOD can be found in the work of Arnold and McLerran ||21| , who studied the 
properties of a dilute gas of sphalerons. They expressed the multiple-sphaleron configura- 
tions as the superposition of many single sphalerons, with partition function approximated 
as above. 

The partition functions in ( [2.10| ) can be evaluated by the saddle-point method, expanding 
the lagrangian field in ( [2.4J) as <f>(x,r) — > <Pb(%) + r l(.%i T ) f° r Z(Pb) an d 4>{^i T ) ~ ¥ Vf + 
((x,t) for Z((pf). r)(x,T) and ((x, r) are small perturbations around the classical field 
configurations ifb{x) and (pf, respectively. These perturbations around each configuration 
bring the temperature corrections to the nucleation barrier into the problem, as we shall see. 
Up to 1-loop order one keeps the quadratic terms in the fluctuations tj(x,t) and ((x,t) in 
the expansion of the scalar field in the partition function (|2.4|) . This way one can write the 
following expressions for Z((pb) and Z(ipf), respectively, 
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l—loop order 



S B (<Pb) 



Dr] exp 



dr / d 3 x-r] [-D E + V"(<p b )] V 



(2.11) 



and 



%>/) 



| DCexpl-^dr / d 3 ^([-n B + y> / )]c] , (2.12) 



where V"(ip) 



d 2 V{4>) 



and D E = J^j + V 2 



Performing the functional Gaussian integrals in (|2.11|) and ( 2.12 ) one gets the following 
expression for the ratio between the partition functions, > a PP earm g m ( P-ICQ : 



Zi^ipb) 1— loop order 

Z((p f ) 



det(-n E +V"(<f b )), 



-AS 



det(-n E + v»&f)h\ " ' (2 ' 13) 

where [det(M) /3 ] _ 5 = J exp { - J"<f / rf 3 x|r/[M]r]} and AS = S E (</? 6 ) - S E (<Pf) is the 
difference between the euclidean actions for the field configurations and Note that 
S E (ip), and hence AS", does not include any temperature corrections. 

From ( p.!0| ) and ( |2.13| ), the free energy of the system, T = — In Z can be written as 

~det(-n E + V"(<p b ))i 



T 



-T 



-AS 



(2.14) 



_det(-D E + V"(tp f ))p_ 
As is well-known, the determinantal prefactor evaluated for the bounce configuration 
has one negative eigenvalue, signalling the presence of a metastable state, and also three 
zero eigenvalues related with the translational invariance of the bubble in three-dimensional 
space. Because of the negative eigenvalue, the free energy T is imaginary. As shown by 
Langer |J (see also Ref. fLlf), the decay rate is proportional to the imaginary part of T 



K 



TFT 



-ImjF 



(2.15) 



where \E_ | is the single negative eigenvalue. In general it depends on non-equilibrium aspects 
of the dynamics, such as the coupling strength to the thermal bath. 



III. EVALUATION OF THE DETERMINANTS 



In this Section we compute the ratio of the determinants appearing in the decay rate, 
and show how it provides a finite temperature correction to the nucleation barrier. 
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First recall that for static field configurations AS is given by 

AE 



//\ H 
d 3 xi£ E (v b )-£ E (v f )] = — , (3.1) 

where AE is simply the nucleation energy barrier, that is, the energy of a critical nucleation 
within the metastable phase. For example, in the thin- wall approximation of Eq. (|2.8|), AE 



is 

AirR 3 

AE = — AV + A7iR 2 a , (3.2) 

o 

where o"o is the tree-level surface tension of the bubble wall (i.e. with no corrections due to 
fluctuations around the bubble wall field configuration <f wa ii) 

r+AR 

ct - / dr[C E ((p wall ) - C E ((pf)} ■ (3.3) 

J-AR 

Using (|3.1|) in Eqs. ( [2.14|) and (|2.15| ) we can write the decay rate as 



1 ( 


AE\ 


exp 1 


T ) 



det(-a E + V"( Vf )),_ — » - ' ■ (3 ' 4) 
The determinantal prefactor in ( |3.4| ) will provide the temperature corrections to the nu- 
cleation barrier. This should come as no surprise, given that the determinant is obtained 
from integrating over thermally induced fluctuations about the tree-level bubble configura- 
tion. (Recall that we are only considering the regime in which thermal fluctuations are much 
larger than quantum fluctuations.) Once the negative and zero eigenvalues are taken care of, 
the positive eigenvalues are easily associated with entropic contributions to the activation 
energy due to stable deformations of the bubble's shape, as in classical nucleation theory. 

We now proceed to show how to incorporate temperature corrections to the nucleation 
barrier. This can be done without an explicit evaluation of the eigenvalues of the operators in 
the determinants. As we show next, all that we need is to separate consistently the positive 
eigenvalues from the negative and zero eigenvalues, and then show how the former can be 
exponentiated. In principle, the computation of the determinantal prefactor in ( |2.14|) can 
be performed by two different methods. The first involves obtaining directly (analytically, 
or more realistically, numerically) the eigenvalues for the determinants in Eq. (|2.14|) . The 
second method consists in developing a consistent perturbative expansion for the ratio of 
the determinants. We now examine both these possibilities. 
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A. Evaluating the Determinantal Prefactor: Eigenvalue Equations 



Consider the eigenvalue equations for the differential operators that appear in the deter- 
minantal prefactor, 



[-n E + V(<p f ))* f (i) = e%i)<s> f (i) 



and 



(3.5) 



-n E + v"(tp b )]$ b (i) = sl(i)® b (i) . 



(3.6) 



In momentum space one writes, e 2 = u 2 + E 2 , where u n = 2j ^ L , n — 0, ±1, ±2, . . ., for 
bosons (for fermion fields u n = ( 2n + 1 ) 7r ) - From (|3.5|) and ( p.6|) one can write the determinant 
ratio in Q2.14 ) as 



det(-n E + V"(y f ))t 
det(-n E + V>'(tp b ))c 



exp I - In 



exp { - In 



det(-n E + V»(<p f ))i 
det{-n E + V"(<p b ))i 

n^n, (u 2 n + E 2 \ 
n+r-ocn.-K+^O')) 



(3.7) 



Using the identity, 



n 1 + 3 



sinh(7rz) 



n=l 



7T,2 



(3.8) 



and taking into account that we have in the denominator of ( ^.7[ ) one negative and three 
zero eigenvalues, one obtains (for details see Appendix A) 



det(-n E + V"& f )) f 
det(-n E + V"(<p b )) 



V 



T 1 



z|£_| sin (^) 

p 



AE 



2txT 



x exp 

- E ' 



E b (j) + In (l - e 



(3.9) 



The factor V 



AE 



2ttT 



in the right hand side of Eq. (|3.9| ) comes from the contribution of the 
zero eigenvalues, which can be handled as in ref. ||, through the use of collective coordinates 
corresponding to the position of the bubble. V is the volume of three space, while the prime 
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in J2j is a reminder that we have excluded the negative and the three zero eigenvalues from 



the sum. Thus, the argument of the exponential incorporates only the contributions from 
the stable vibrational modes on the bubble, very much like in Langer's result |J. 

Substituting Eq. ( |3.9|) in (|3.4j) one obtains the following expression for the nucleation 
rate per unit volume, T = 71/ V, as defined in Eq. Q2.15Q : 

AF(T)~ 



T = „4T 4 exp 



where we have denoted by A the dimensionless factor 



(3.10) 



A 



1 



\E-\ 
2T 



71 



sin (J^l) 



AE 



2nT 



(3.11) 



Note that the zero-mode contribution in the prefactor depends on the energy barrier of the 
critical nucleation, and not on its free energy as in Eq. (El). In Eq. (ft.lC ) we have incorpo- 
rated the exponential contribution from the prefactor into the definition of the temperature 
corrected nucleation barrier, which we call AF(T). Within the thin-wall approximation we 
can write 

AttR 3 



AF(T) 



-AV cS (T)+4nR 2 a(T), 



(3.12) 



where 



and 



AV e s(T) — V(<ff) — V((f t ) + T 



d 3 k 
(2tt) 3 



In 



-f}yjk 2 +-m?{ip f ) 



- T 



a(T) = a + 



E 



d 3 k 
(2^)3 



T 



In 



AttR 2 


2 



^E wall (j)+ln(l-e-^ 



u(j) 



§%(i)+ln(l-e- 



(3.13) 



(3.14) 



In Eq. ( |3.13| ) we have substituted the discrete sums by integrals over momenta. [For field 



configurations tpf and (pt, we have the continuum eigenvalues, E'j = k 2 + m 2 (ipf) and E 2 



k 2 + m 2 ((pt), respectively, with m 2 ((pf) 



d 2 VW 



=<Pf 



and m 2 (ip t 



d 2 V(4>) 



=<pt-l 



We have 



12 



omitted the usual zero-temperature ultraviolet divergent terms, / d 3 kyk 2 + m 2 (</>), from Eq. 
( p,13| ) since they can always be handled by the introduction of the usual counterterms |22| . 

In Eq. ( p. 14 ) E wa u(j) are the eigenvalues related with the bubble wall field configuration 
(fwaii- Thus, within the thin- wall approximation, the problem is reduced to the computation 
of these eigenvalues for a field configuration describing the bubble wall, a non-trivial task. It 
is possible, however, that the temperature corrections to the surface density are negligible. 
For example, in the context of the QCD transition, the surface density is obtained from 
lattice calculations, and is shown not to be very sensitive to temperature |~2] . 



It is easy to see from Eq. ( |3.13| ) that AV c q(T) is the usual 1-loop approximation to the 
finite temperature false vacuum energy density [^2] . The second term in the right hand side 
of (3.14 ) comes from the finite temperature 1-loop contribution to the surface tension cq, due 
to thermal fluctuations on the bubble wall. Thus, by exponentiating the contribution from 
the stable modes, the activation energy for the critical bubble becomes indeed an activation 
free energy. Note however, that contrary to Eq. ( |1 . 1| ) , the free-energy functional is evaluated 
for the tree-level bounce. In Section 4 we will compare the results obtained with the two 
approaches. 



B. Evaluating the Determinantal Prefactor: Perturbative Expansion 



A second approach to the computation of the determinantal prefactor in ( |3.4|) consists 
in developing a perturbative expansion for it. First we write the determinantal ratio as 



det(-n E + V"(ip f ))p 
det'(-a E + V"(<p b )) p 



= exp|-Tr In [-D E + V"(tp f )] p - 
~ ^Tr' \n[-n E + V" , 



(3.15) 



where we have used in ( |3.15| ) the identity In det M = Tr In M and the prime in both sides 
denote that the negative and the zero modes have been omitted. (They are treated as in 
previous Section.) 

We rewrite ( |3.15|) as 
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det(-D g + V"(y f ))p 
det'(-D E + V"(<p b ))f, 



cxp <j — -Tr In 



1 + Gp( Vf ) \V\ Vb ) - V"{ Vf )\ 



where 



-D E + m 2 ((p f ) 

is just the propagator for the scalar field 0, with m 2 (ipf) = V"(ipf). 
Expanding the logarithm in Q3.16 ) 



m 



powers 



Gp(ipf) \V"((pb) — V'i^Pf)], we obtain the graphic representation, 



(3.16) 



(3.17) 



of 



Tr \n{l + G p { Vf )[V"{ Vh )-V"{ Vf )]} 






+ 



(3.18) 



where the dashed lines correspond to the background "field" [V"(ipb) — Viff)} an d the 
internal lines denote the propagator Gp{(pf). The expression ( |3.18| ) can be written as 

+oo j i\m+l 



Tr ln{l + Gp((ff) [V"(cp b ) - V"(cp f ))} = J2 



m=l 

+oo 

x E 



m 



d*x[V"{y b )-V"{y f )) m x 



d 3 k 



(27T) 



The sum in m can be performed and one obtains 
Tr ln{l + Gpiiff) [V"M - V"(ip f )}} = / d 3 x £ 



d 3 k 



(2vr 



u 2 + k 2 + m 2 (iff) 

ln l , v"M - v"(<P f ) 



(3.19) 



uj 2 + k 2 + m 2 (ipf)_ 

(3.20) 



This expression can be further simplified by means of the identity, (E 2 (k) = k + m 2 (ipb), 



and E 2 f (k) = k + m 2 (ip f )) 



E ^ 



<4 + tf(k) 

co 2 n + E%k) 



rt+21n(l-e-^ 
PEf(lt)-2hi(l-e- f>E f® 



(3.21) 



The terms proportional to (3 can be renormalized. We are left with the familiar tempera- 
ture corrected effective potential (we neglect terms coming from zero temperature quantum 
corrections) |22| . 
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Vrf(<p,T)=V(<p)+T 



d 3 k 



In 



^ _ e -/3y/P+m?{<p) 



(3.22) 



Within the thin-wall limit, using ( |3.21| ) into ( |3.20|) , and substituting into Eq. ( |3.16| ), we 



obtain, from Eq. (3^) the temperature corrected barrier, 

AttR 3 



AF(T) 



-AV cS (T)+4irR 2 a(T) 



where 



AV eS (T) = V(( Pf )-V(( P t)- 



d 3 k 

(3 J (27P 



In 



-PE f {k) 



(3.23) 



(3.24) 



and 



1 r 1 + °° r 



d 3 k 



In 



1 + 



oj\ + k 2 + m 2 (ip f ) 



(3.25) 



2(3^ J (2vr)3 

Expression (|3.24f ) for the temperature corrected false vacuum energy density can be easily 
identified with AV e s(T), as defined in Eq. ( |3.13| ). Also, expression ( |3.25| ) for cr(T) can be 
identified with the temperature corrected surface tension of the bubble wall, by writing it 

as 



a(T) 



AnR 2 



[F{y walh T)-F{y h T)] 



with 



F(^,T) = -ilnZ(^) 



(3.26) 



(3.27) 



where F(<p,T) is the free energy for the field configuration (p. cr(T), as given by ( |3.26| ), is 
then the free energy difference (per unit area of the bubble) between the two configurations 
(Vwaii and <ff), which defines the surface tension. As shown in Ref. |23], the expressions for 
the free energies in (|3.26|) , at 1-loop order give ( p.25| ). 



IV. EFFECTS OF INTERACTIONS ON THE DECAY RATE 

In the previous two Sections we have obtained the temperature effects on the nucleation 
barrier by assuming that the system becomes trapped in a metastable phase as the temper- 
ature drops below the critical temperature Tq. This assumption is not very realistic, and 
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was addopted so that we could stress the main points of the calculation without having to 
worry about the effects of interactions of the order parameter with other fields. However, 
in models of interest within a cosmological context, such as the standard electroweak model 
and some of its extensions, the Universe becomes trapped in a metastable phase due to the 
interactions of the Higgs field (or effective scalar order parameter) with other massless fields, 
such as gauge bosons or fermions. In fact, it is possible to formulate this statement in terms 
of conditions that a general effective potential must satisfy, if it is to develop a metastable 



phase below Tq |24| . Basically, the condition states that the mass gap between the symmet- 
ric and broken-symmetric phases must be large enough so that massive fluctuations away 
from the symmetric minimum are suppressed. Since the mass gap is given in terms of the 
(temperature dependent) vacuum expectation value of the Higgs field and of its coupling to 
other fields, the condition states that the transition cannot be too weakly first order. 

The question then is what effective potential should be considered when calculating the 
decay rate. As recently shown by E. Weinberg for the zero-temperature Coleman- Weinberg 
models (which exhibit symmetry breaking only due to radiative corrections), the effective 
potential relevant in the calculation of the bounce solution is obtained by integrating over 
all fields but the scalar field; the radiative corrections coming form these fields modify the 
vacuum structure of the model, making metastability possible. In this Section, we argue 
that the same procedure must be followed when calculating the nucleation rate at finite 
temperature. We first review, in general terms, how to take into account the interactions of 
the scalar field with other fields in the calculation of the decay rate. (For details see, Ref. 
T8fl.) Then we apply our method to an example involving fermions, explicitly comparing 



our results to those obtained from Eq. (|1.1|) . 

A. General Formalism 

Let us consider a system described by a scalar field and a set of fields £ (bosonic or 
fermionic fields) which are coupled to <fi. The partition function of the system is given by 

Z = J D(j)Die- SE( - 4 '' i) , (4.1) 
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where Se(4>,£) is the euclidean action of the system and the functional integration is carried 
over field configurations subject to periodic boundary conditions, (f)(x, 0) = 0(x, (3), for 
bosons or antiperiodic, i/j(x,0) = —ip(x,f3), for fermion fields. 

If one integrates out the £ fields in (|4.1| ), the partition function can be written as 



fe- wB W , (4.2) 
where 

W p {(j)) = —In / Dfe -5 *^' . (4.3) 
W^(4>) can be viewed as an effective action for the scalar field 0, where only £-loop terms 



are included. Note that these £-loop terms introduce finite temperature corrections in (4.3). 



Se(4>,£) in ( |4.1|) and (|4.3| ) includes renormalization counter-terms and, if one or more of the 
^-fields is a gauge field, Se{(^, £) also includes the gauge fixings and the corresponding ghost 
terms. 

The procedure is now, in principle, straightforward. We evaluate the partition function 
in Eq. Q4.2| ) semi classically by expanding the effective action W^{4>) around its extremum 
configuration, which will be the bounce. Note that the bounce will include the radiative 
corrections coming from the fields that couple to <ft, but not from <fi itself. The determinan- 
tal prefactor can be evaluated as before, by considering the negative and zero eigenvalues 
separately from the positive eigenvalues. However, as pointed out by Weinberg ||18|| , W^{(f>) 
cannot always be obtained in closed form, being in general a nonlocal functional. He pro- 
poses to resolve this difficulty by considering a local action Wo((p) which is close enough to 
the original one. From a derivative expansion of W^{4>), Wo((p) is found to be 

W (<f>) = J d A x E [i(<9^) 2 + Vl_ ioop (0)] , (4.4) 

where Vi-i oop ((f)) includes the tree level potential V((p) and the 1-loop contributions coming 
only from the £-field integration. The bounce solution can then be obtained from ( [4.4|) . 
Wo(4>) should be a good approximation to W^(4>) as long as the typical interaction length 
scale set by the field(s) £ is shorter than the scale for the nonlocality of W 13 ((f)). 
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In what follows we will apply the above formalism to a specific example of a scalar field 
coupled to fermions. For large enough Yukawa couplings, this model has been shown to 
satisfy the conditions for metastability specified in Ref. |24|. 

B. Application: Real scalar field coupled to fermion fields 

Consider a model of a real scalar field (f) coupled to fermion fields with a lagrangian 
density 

C = l -{d^f - V{4>) - gm W> , (4-5) 



with V{4>) given by 



V(4>) = ^ - f 3 + , (4.6) 



where A and a are positive, dimensionless constants and fi 2 > is a (mass) 2 parameter. 
A large enough coupling to fermions guarantees that a metastable phase is possible as the 
system is cooled below Tc; the high temperature minimum of the 1-loop effective potential, 
(<p)t — x+ig 2 f 1 ' ^ es ^° ^ ne °f ^ ne max i mum °f the potential at T = Tq, and thermal 
fluctuations away from the symmetric minimum are suppressed. Numerical values for the 
couplings satisfying these conditions can be found in Ref. |24|. 
The partition function for this model is given by 

' D(fiDipDipe~ to *rf**mi>,® . (4.7) 

As the fermion fields appear quadratically in ( |4.7|) , they can be integrated out, giving 

Z = J D(j>e- w °W , (4.8) 

where 

W Q {4>) = J^dr J rf 3 x[i(9 r 0) 2 + i(V0) 2 + V(0)] -Trln(- - ig</>) p , (4.9) 

and 
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THn(- p-igflp = lndet(- - ig<j>) . (4.10) 



Thus, Wq{4>) can be written as 



Wo(<p) = dr d A x 



o 



\(d T <p) 2 + l -(y<p) 2 + v^T) 



(4.11) 



where, for sufficiently smooth fields (see Appendix B), the effective potential obtained af- 
ter integrating over the fermions and renormalizing is (we will drop all zero temperature 
quantum corrections) 

^(0,T) = n0)-4T/^ln(l + e^v^W) , (4.12) 

The temperature dependent term accounts for finite temperature corrections coming from 
fermion loops. V^((p,T) is the potential we should use to compute the bounce. Note that, 
neglecting the T = quantum corrections, the high-temperature limit of V^((f), T) is approx- 
imately 

^(0,T)~V(0) + ^ 2 2 , (4.13) 

so that for (T//i) 2 > (9a 2 /X — 24) /4g 2 , a condition which is easily satisfied for reasonable 
values of the couplings, the high-temperature minimum is (v?)t — 0. 

Once we have the action Wq(4>), the bounce is obtained as a solution of 



0. (4.14) 



Thus, for a static, spherically symmetric configuration, the bounce configuration y9&(r) will 
be a solution of 

d 2 4) 2d<p fru.rris fA * K \ 



dr 2 r dr 



with boundary conditions, lim r _ >0O ipb(r) = (pf ~ and ^|r=o — 0. (From now on (pj and 



(ft should be understood as being the minima of ( 4.12|) and not the minima of the tree-level 
potential V((p).) 

The procedure is now identical to that of Sections 2 and 3. Having a bounce solution 
we can evaluate the partition function written in Eq. ( [4.8D semiclassically, exactly as was 
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done in Eqs. ( |2.11|) and ( |2.12| ), by expanding around iff and ip b . We then obtain, from Eqs. 
( p.l4j ) and (|2.15|) , the nucleation rate, 

1 

'det(-n E + m^((p b )) fi ' 



where m 2 Jif) 



K = 



\E. 



-Im 



-AWo 



7T [det(-n E + m 2 p(<pf))p_ 
=(p , with V^(4>) obtained above. AW is given by 



(4.16) 



(4.17) 



where Wq(4>) was defined in Eq. (|4.11|) . 

In order to proceed, we must rewrite the determinantal prefactor explicitly isolating the 
negative and zero modes from the positive modes. This is done following the same steps 
of Section 3, although now we must handle the fermionic contribution to the determinants. 
The details of the perturbative expansion for the fermionic determinantal prefactor are given 
in Appendix B. We can then write the nucleation rate per unit volume as 



T 4 \E-\ 
Y = 221_ 



\aw - 


3 
2 


AF(T)~ 


[ 2tt . 


exp 


T 



(4.18) 



* sin (Jf£l) 

where AF(T), the bubble activation free energy in the 1-loop approximation, is given by 



1 1 +0 ° 

AF(T) = I d 3 x{- (Vipb) 2 + - MVf) + 2p E 



As usual, the sum over n can be performed and we get, 

1 



d 3 k 



In 



ul + k 2 + m 2 fi (<p f )_ 
(4.19) 



AF(T) 



d 6 x 



(Vwy + V eS (<p b ,T) - V eS (<p f ,T) 



(4.20) 



where the effective potential Ves((f),T), is given by (neglecting zero temperature quantum 
corrections) 

d 3 k 



V eB (<f>,T) = V(<f>)-4T 
d 3 k 



+ T 



In 



In f 1 + e -/JV^W 
(2tt) 3 V 

1 _ ^-l3^P+mlw\ 



+ 



(4.21) 
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where V(4>) is the tree level potential ( |4.6|) and the mass term appearing in the scalar loop 
contribution is ~ V"((f)) + ^-g 2 in leading order in the fermion loops. It is instructive 

to contrast this result with that obtained for self-coupled scalars, Eq. (|3.22| ). The coupling 
to fermions modifies the scalar mass propagating in the loops. This effect naturally improves 
the infrared behavior of the theory, and can be of importance in weak first-order transitions. 
We will say more about this later. Note also a crucial difference between this expression 
for the free energy barrier and the expression for the free energy barrier in Eq. fll.lj) : Here, 
the bounce is obtained with the effective potential that does not include the corrections 
coming from scalar loops. The corrections from scalar loops which appear in the last term 
of Eq. (14.21 ) are thermally induced fluctuations about the bounce solution computed with 
V^/,(</>, T). In the usual expression for the nucleation barrier the bounce is obtained from the 
full effective potential including the scalar loops. The two expressions are definetely not 
equivalent, even though for small scalar self-couplings the differences are negligible. In order 
to illustrate the differences let us look at a specific example. 

In Fig. 2 we contrast the two approaches by comparing the nucleation barriers as a 
function of the temperature for a fixed set of coupling constants. The barriers in the figure 
were obtained by a numerical integration of the bounce equation including the relevant loops 
according to each approach. For clarity let us call the nucleation barrier obtained in the usual 
approach, i.e., by including the scalar loops in the bounce calculation, the scalar barrier. 
The nucleation barrier obtained without including the scalar loops we call the fermionic 
barrier. In Fig. 2 we take A = 1.0, a = 2.0, and g 2 = 0.5. Since A controls the strength of 
the scalar corrections, we expect the differences between the two barriers to be noticeable. 
We find that this is indeed the case, noting that as we approach the critical temperature 
(that is, as we move closer to the thin-wall limit) the differences between the two barriers 
increase, with the scalar barrier always larger than the fermionic barrier. This is precisely 
what one expects if the scalar corrections are entropic corrections to the nucleation barrier. 
Thus, the nucleation barrier used in expression Eq. (|1.1|) is overestimated for large enough 
scalar corrections. 
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Finally, we point out two additional differences between the results. First note that the 
contribution from the zero modes to the prefactor depends on AW o (0), as opposed to AF(T). 
This could be important for weak transitions in which the prefactor may play a relevant role. 
Most importantly, the expression for V e e((f>, T), Eq. (|4.21|) , differs from the usual 1-loop finite 
temperature effective potential by the mass term for the scalar field loops, rrii((j)) = V^(4>). 
Since we have used the stationary points of W ((f>), Eq. ( |4.11| ), as opposed to the stationary 
points of Se(4>, ij))-, as the effective "background" fields in the saddle-point evaluation of the 
partition function, the scalar field propagator carries the finite temperature mass m^(0). 
The propagator is dressed by the quantum corrections due to fermion loops. In the usual 
1-loop finite temperature effective potential, the stationary points are obtained from the 
tree level action, with mass term for scalar loops, 77?.q(0) = V"(4>). This results in the usual 
negative mass terms related to the change in convexity of the effective potential between 
the inflection points, and, in the case of very shallow potentials, in bad infrared behavior 
near iff. The incorporation of the fermionic corrections to the scalar propagator, which is 
demanded by our method of calculation atenuates these problems. In the example above, 
the scalar mass gets dressed by fermionic loops, being given by m 2 n{4>) ~ V"((f>) + ^g 2 T 2 , 
where V(<fr) is the tree-level potential ( |4.6| ). The temperature term in mp works as the 
infrared regulator for small values of m^cp) = V"(<ft). This result is independent of the 
particular model studied. Similar conclusions have been obtained in Ref. [13| for scalar 
electrodynamics. 



V. CONCLUSIONS 

In this paper we examined in some detail the computation of false vacuum decay rates 
at finite temperatures in the regime in which quantum fluctuations are negligibly small 
compared to thermal fluctuations. We have shown that temperature corrections to the 
nucleation barrier can be obtained from a saddle-point evaluation of the partition function 
in a dilute gas approximation. In fact, the temperature corrections are simply due to the 
positive eigenvalues from stable fluctuations around the critical bubble. That is, they are 
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the entropic contributions due to thermally induced deformations on the bubble. 

Even though this result has been known in classical statistical mechanics for more than 
two decades ||, we believe that a consistent treatment within field theory is still lacking. 
Although we left many questions unanswered, we hope to have clarified some of the issues 
involved in the calculation of finite-temperature decay rates. Of particular importance is the 
fact that the bounce is not obtained from the full 1-loop corrected effective potential, but 
from the potential excluding the scalar loops. Thus, for a self-interacting scalar, the bounce 
is obtained from the tree-level potential. The full finite temperature potential appears in the 
exponent only after properly accounting for the positive eigenvalues of the determinantal 
prefactor. That is, the scalar contributions account for entropic corrections to the nucleation 
barrier. We obtained a temperature corrected nucleation barrier which can differ from the 
usual result. We showed this to be particularly true for sufficiently large scalar self-couplings 
in the vicinity of the critical temperature for the transition. 

Also, we found that the interaction with other fields gives rise to a potential which is 



better behaved in the infrared. (See also Ref. ||13|| .) This result is the finite-temperature 
equivalent to what E. Weinberg found for the zero-temperature case, once the integration 
over the other fields is performed [|1^ . 

The reader may be wondering if our results will have any consequences to current work on 
the electroweak phase transition. The answer depends on the Higgs mass. For a sufficiently 
light Higgs it is consistent to neglect the contribution from scalar loops to the effective 
potential. In this case, the usual estimate for the nucleation barrier is a valid approximation. 
However, the situation may change for a heavier Higgs. Given that the experimental lower 
bound on the Higgs mass is now above 60 GeV, we believe it worthwhile to study this 
question in more detail, keeping in mind that the transition becomes weaker as the Higgs 
mass increases. 
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Figure Captions 



Figure 1: A typical asymmetric double- well potential. 

Figure 2: A comparison of the nucleation barrier as a function of temperature, in 
units of mass parameter /i, obtained by including (stars) and excluding (dots) scalar 
loops in the computation of the bounce. The parameters in the tree-level potential are 
A = 1.0, a = 2.0, g 2 = 0.5. 

Figure 3: A comparison of the terms g 2 (j) 2 (r) and g d(p(r)/dr appearing in Appendix B. 
The parameters in the tree-level potential are A = 1, a = 0.56, g 2 = 0.5. 



APPENDIX A: 



Let us show that the determinant ratio of Eq. ( |3.7| ), denoted here as R, gives (|3.9|). 
Separating the negative and zero eingevalues in the denominator of Eq. ( |3.7| ), one can write 

n; 



R 



cxp { -In 



-oo 
in=— oo 



Yli(ul + E 2 (i\ 



it 



■l + El){ul + EiytfAul + El(j)) 



(Al) 
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where the prime in rL means that the negative eigenvalue, E^_, and the three zero eigen- 
values, Eq, are now excluded from the product. The term for n = in {uj\ + Eq), can be 

3 

handled as in Ref. 0, resulting in the factor V |MI 2 in Eq. (|3.9|) . Separating the n = 



modes both in the numerator and the denominator of (|AT|), and using the identity (|3.S| ), we 

get, 



R = V 



AE 



exp 



4 + E-E']lnn^-ln(^) 1/2 -ln 



2vrT 

I E ' - E I h) ; E 



n=l 



sin(§|£_|) 



2 \- LJ - 



+ 



§%(0+ln(l-e-^« 



+ 



E' 



^ 6 (j) + In (l - e H»*« 



(A2) 



In the above expression we used that the negative eigenvalue can be written as (£^.) 2 
il-E 1 -!- Remembering that Yfj has four eigenvalues less than we can write 

,j/2 , "sin(f \E_\y 



R = V 

+ E 



AE' 


1 ( 
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2nT_ 


exp < 



+ 



|%(i)+ln(l-e-^'« 
which reduces to Eq. ( p.9|) . 



E' 



fs 6 0*)+ In (l-e- 



■0E b (j) 



(A3) 



APPENDIX B: 



In ( [4. 16 ), the exponential term AW$ can be written as (from ( 4.17|) and ( [4.9| ) 



AW = [S E (<p b ) - Trln(- - ig Vh ) p ] - [S E (<p f ) - Tr ln(- 



where = J? dr J d 3 x \Ud^) 2 + V(</>) 



(Bl) 



is the classical action for the scalar field as in 



3|). The exponential in ( |4.16| ) can be written then as 

e -Aw = [ det(- 

|_det(- @-ig<Pf)p\ 



-AS 



(B2) 



with A£ = S^(</? 6 ) - S E (ip f ). 
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If one uses the identity (which follows from charge- conjugation invariance): 



[det(— ft — ig(p)] 2 — det(— ft — ig<p). det(— ft + ig<p) 



det 



D E + g ip )l 4x4 - ig^d^ip 



(B3) 



where 1 4X 4 is the 4x4 unit matrix, then, for ip = ipf, one gets 



det(- ft - igVf)p = det(-D B + g <p f ) p 



(B4) 



For the determinant involving the spherically symmetric bounce ipb, the Dirac matrix j E in 
is radial and one can write 

/ 



(B5) 



V 



1-2x2 
— 12x2 

where 1 2X 2 denotes a 2 x 2 unit matrix. Then, for <p b = <pb( r ), det(— ft — ig(p b (r))/3 can be 
written as 



det(- ft-igy b {r)) p = det tt i+) {<p b ). det Q { -\ip b ) 



(B6) 



where 



& ±) {<p b )=-n E +g 2 ( pl±g 



chph 
Or 



(B7) 



Therefore, the determinants in (|B^ ) can be written as (using that In det M 
TrlnM ) 

det(- ft- ig<pb)/3 



det(- ft-igtpf) 



exp {Tr ln(- ft - igip b )p - Tr ln(- ft - igpf)p} 



exp I Tr In 
2Trln 



■^E + gVb + 9^ 



+ Tr In 



-1/3 
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E + g 2 if 2 f 



-1/3 

(B8) 



As in Q3.160 , one can write ( |B8| ) as 
det(- ft-ig(pb)f3 



det(- ft - igcpf) 



exp < Tr In 



1 + S p {ipf) 



9 &b -<Pf) + g 



di£b 
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+ Trln 



1 + S fi {<p f ) 
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(B9) 
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where, in analogy to Eq. ( |3.16| ), we introduce the propagator 



(BIO) 



-U E + g*v) 

The argument of the exponent in the rhs of (|B9| ) can be written as a series, analogously to 
Eq. flOl): 



Tr In 
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where uj t . 



(2n+l)w 



. As before, ( |B11| ) can be expressed as a graphic expansion similar to 



( |3.1S| ), with the propagators Gp(<pf) replaced now by Sp(cpf) and the external lines given by 
9\ti ~ ¥>}) + 9^ or g\rt - V }) - . 

The determinant factor in (|4.16|) , coming from the functional integration of the scalar 
field, can be evaluated by the same methods of Sec. 3. In (|4.16|) , the determinant term 
det[— O e + mp((pb)]p, with 7n|(<^&) = V^'(</?&), has a negative eigenvalue, E^_, associated 
with the instability of the critical bubble, and the three zero eigenvalues, associated with 
the translational invariance of the bubble. These eigenvalues can be handled as usual, 
giving the preexponential term in ( [4.18| ). The part of the determinant involving the positive 
eigenvalues can be written as an expansion exactly as in ( |3.19| ), 



' det'(-D g + V$(<p b ))p 
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and 
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{^f[ul + P + ml( V} )] 

The sum in m in both ( B 1 1|) and (|B13|) can be performed as in Eq. (|3.19| )). Therefore, from 
Eqs. ( P2|) , (|B11|) and ( B13|) , we can write the relevant part of Eq. (|4.16| ) as 



27 



det(-n E +V»(p b )) 
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(B15) 



+ P + 

Apart from the derivative terms dip^/dr, the momentum integral reproduces the finite 
temperature corrections to the the tree-level potential appearing in AS*. When we wrote 
the expression for AF(T) in Eq. ( 4.19Q , these terms were not included in the effective 
potential Kg (0, T). There are two reasons for negleting this term. First, due to the graphic 
expansion we used for the determinants, it is easy to see that at least at the tadpole level, 
their contribution cancels. Since the tadpole gives the dominant temperature contribution 
to the potential, terms that depend on dtpb/dr will be sub-dominant. Second, it is possible 
to explicitly compare the terms g 2 ^\ and gd(fi,/dr, by obtaining (pb(r) numerically. We 
have performed this comparison for the same set of parameters used in Figs. 2 and 3, and 
convinced ourselves that the derivative term will indeed be sub-dominant. A typical example 
is shown in Fig. 4. Thus, neglecting the term d(fb/dr, we can use Eqs. ( pi4|) and (|B15| ) to 
obtain the expression for AF(T) in ( 4.19Q . 
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